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The  Mediterranean  Sea  produces  a  salty,  dense  outflow  that  is  strongly  modified  by  en¬ 
trainment  as  it  first  begins  to  descend  the  continental  slope  in  the  eastern  Gulf  of  Cadiz. 
The  current  accelerates  to  1 .3  meters  per  second,  which  raises  the  internal  Froude  number 
above  1 ,  and  is  intensely  turbulent  through  its  full  thickness.  The  outflow  loses  about  half 
of  its  density  anomaly  and  roughly  doubles  its  volume  transport  as  it  entrains  less  saline 
North  Atlantic  Central  water.  Within  100  kilometers  downstream,  the  current  is  turned  by 
the  Coriolis  force  until  it  flows  nearly  parallel  to  topography  in  a  damped  geostrophic 
balance.  The  mixed  Mediterranean  outflow  continues  westward,  slowly  descending  the 
continental  slope  until  it  becomes  neutrally  buoyant  in  the  thermodine  where  it  becomes 
an  important  water  mass. 


The  deep  ocean  is  filled  with  water  masses 
that  originate  in  marginal  seas  where  a  cold 
or  dry  continental  climate  can  produce 
dense  water  by  intense  and  sustained  air-sea 
interaction  (1).  This  dense  water  forms  bot¬ 
tom  currents  that  flow  out  from  the  marginal 
seas  and  settle  in  the  open  ocean  at  levels 
that  are  determined  by  their  relative  densi¬ 
ties  (2,  3).  The  temperature,  salinity,  and 
other  chemical  properties  of  these  outflows 
set  the  climate  of  the  deep  sea,  which  is 
known  to  have  been  significantly  different  in 
other  epochs  (4)  from  what  it  is  today.  One 
step  toward  better  interpretation  of  the 
record  of  past  deep-sea  climate  change  or 
toward  prediction  of  future  climate  change  is 
to  understand  what  processes  set  the  proper¬ 
ties  of  these  marginal  sea  outflows. 

We  examined  the  outflow  that  begins  in 
the  Mediterranean  basin  where  the  loss  of 
freshwater  by  evaporation  causes  Mediterra¬ 
nean  seawater  to  have  a  high  salinity,  S  = 
38.4  psu  (practical  salinity  units,  nearly 
equal  to  parts  per  thousand),  and  conse¬ 
quently  a  high  density,  <x9  =  28.95  kg  m~3 
(5).  This  density  exceeds  that  of  North 
Atlantic  bottom  water,  cre  =*  27-95  kg  m-3, 
and  is  much  greater  than  the  density  of 
North  Atlantic  surface  water  found  just  out¬ 
side  the  Mediterranean.  The  resulting  hori¬ 
zontal  contrast  in  density  drives  an  inverse 
estuarine  circulation  through  the  shallow 
and  narrow  Strait  of  Gibraltar  (Fig.  1). 
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Fresher  and  lighter  Atlantic  surface  water 
forms  the  inflowing  upper  layer  that  is  trans¬ 
formed  by  freshwater  loss  into  saline,  dense 
Mediterranean  water  that  forms  the  outflow¬ 
ing  lower  layer  (6) .  The  Mediterranean  out¬ 


flow  spills  over  the  shallow  sills  in  the  Strait 
of  Gibraltar  and  cascades  down  the  conti¬ 
nental  shelf  and  slope  into  the  Gulf  of  Cadiz 
and  the  open  North  Atlantic  (Fig.  2). 

Despite  its  initially  very  high  density, 
the  Mediterranean  outflow  water  does  not 
reach  the  bottom  of  the  North  Atlantic 
because  it  entrains  (7)  a  substantial  volume 
of  the  overlying  North  Atlantic  Central 
water  while  still  in  the  Gulf  of  Cadiz. 
Entrainment  reduces  the  high  salinity  and 
density  of  the  outflow,  and  the  resultant 
mixed  Mediterranean  water  that  finally 
reaches  the  open  North  Atlantic  has  a 
density  cre  =  27.6  ±  0.1  kg  m-3  that  makes 
it  neutrally  buoyant  at  depths  of  about  800 
to  1300  m  in  the  main  thermocline  (8-/3). 
Mixed  Mediterranean  water  is  considerably 


Fig.  1.  (A)  North  Atlantic 
salinity  distribution  (in 
practical  salinity  units,  H 
and  L  mark  highs  and 
lows)  on  a  density  surface 
that  is  near  the  salinity 
maximum  of  the  Mediter¬ 
ranean  outflow  (its  depth 
is  1200  m  in  the  tropical 
and  subtropical  North  At¬ 
lantic,  and  it  reaches  the 
sea  surface  in  the  Norwe¬ 
gian  Sea  and  in  shelf  wa¬ 
ters  northwest  of  the  Gulf 
Stream).  A  high-salinity 
tongue  of  mixed  Mediter¬ 
ranean  water  originates 
from  the  northern  Gulf  of 
Cadiz  and  extends  west 
across  the  North  Atlantic 
basin  and  north  toward 
the  Norwegian  Sea.  The 
other  water  masses  at  this 
density  level  all  have  low¬ 
er  salinity  (Labrador  Sea 
water  in  the  northwest, 
surface  water  of  Arctic  or¬ 
igin  in  the  northwestern 
Norwegian  Sea,  and  Ant¬ 
arctic  Intermediate  water 
to  the  south).  [Adapted 
from  figure  2  of  ( 74)]  (B) 
Salinity  in  the  eastern 
North  Atlantic  and  Gulf  of 
Cadiz  at  a  depth  of  1 200 
m.  The  dashed  line  run- 
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ning  through  the  Strait  of  Gibraltar  denotes  the  axis  of  the  outflow  that  is  shown  in  vertical  section  in 
Fig.  2A.  [Adapted  from  (3fi)J  (C)  The  maximum  observed  velocity  of  outflow  currents  in  the  eastern 
Gulf  of  Cadiz  as  measured  by  expendable  current  profilers  (XCPs).  The  lines  of  stations  that  make  up 
sections  are  labeled  A  through  F  in  order  of  increasing  distance  away  from  the  Strait  of  Gibraltar.  Note 
that  the  outflow  current  crossed  isobaths  (dashed  lines,  depth  in  meters)  toward  deeper  water  in  the 
vicinity  of  section  C  and  was  more  or  less  parallel  to  the  local  isobaths  near  sections  E  and  F. 
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saltier  and  warmer  than  thermocline  water 
masses  that  have  subpolar  origins  and  makes 
a  conspicuous  anomaly  that  extends  south¬ 
west  across  the  subtropical  North  Atlantic 
and  north  along  the  eastern  boundary  of  the 
northern  North  Atlantic  (Fig.  1A)  (14). 

A  Field  Survey  of  the 
Mediterranean  Outflow 

We  conducted  an  intensive  field  survey  of 
the  northern  Gulf  of  Cadiz  in  the  fall  of  1988 
to  characterize  the  physical  dynamics  of  the 
Mediterranean  outflow  and,  more  specifical¬ 
ly,  to  learn  where  and  how  the  outflow  is 
modified  by  mixing.  The  field  data  include 
99  full-depth  profiles  of  temperature  and 
salinity,  T(z )  and  S(z),  that  were  made  with 
a  conductivity,  temperature,  and  depth  in¬ 
strument  (CTD)  (15).  At  56  of  the  CTD 
stations,  we  obtained  a  profile  of  the  hori¬ 
zontal  current,  U(z),  with  an  expendable 
current  profiler  (XCP)  (16)  (Fig.  3)  and,  at 
29  of  those  stations,  a  profile  of  turbulence 
intensity  with  an  expendable  dissipation 
profiler  (XDP)  (17)  (Fig.  4). 

About  two  thirds  of  the  CTD  stations 
were  laid  out  in  the  Gulf  of  Cadiz  in  eight 
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Fig.  2.  (A)  A  salinity  section  along  the  axis  of 
the  outflow  made  with  data  from  conductivity, 
temperature,  and  depth  instruments  (CTDs) 
(location  of  this  section  is  shown  as  the  dashed 
line  in  Fig.  1 B).  The  origin  (0  km  downstream)  is 
the  main  sill  at  Camarinal  (5°46'W).  This  section 
is  very  similar  to  that  produced  by  Heezen  and 
Johnson  (8)  with  data  taken  until  about  1965. 
(B)  Traditional  core  value,  S^.  the  maximum 
salinity  within  a  section;  the  transport-weighted 
salinity  within  the  outflow,  S^,  which  we  be¬ 
lieve  is  a  better  measure  of  the  average  salinity 
of  the  outflow;  and  S*,,,  the  salinity  of  the  North 
Atlantic  water  just  above  the  outflow.  The  val¬ 
ues  of  Soy,  and  S*,,  are  nearly  identical  to  those 
shown  in  figure  12  of  (8). 


sections  that  were  oriented  across  the  path 
of  the  outflow  (labeled  A  through  H  in  Figs. 
1C  and  2A).  We  found  these  data  to  be 
similar  to  hydrographic  data  acquired  in 
similar  surveys  that  have  been  conducted 
since  the  1960s  (18). 

We  estimated  the  transport  of  volume, 
salinity,  and  density  through  each  of  the 
sections  (Fig.  5)  by  combining  the  CTD 
and  XCP  current  data.  By  analyzing  the 
transport  data,  we  calculated  a  transport- 
weighted  average  salinity  for  each  section, 
STwtd,  defined  as  the  salinity  transport  in 
the  outflow  divided  by  the  volume  trans¬ 
port,  which  should  give  a  clear  and  repre¬ 
sentative  measure  of  salinity  variation  along 
the  outflow  path  (Fig.  2B).  As  expected, 
STwtd  has  an  initial  value  of  38.4  psu  in  the 
Mediterranean  Sea  and  decreases  gradually 
as  the  outflow  crosses  the  main  sill  at 
Camarinal  (5°46'W).  A  corresponding 
gradual  change  is  evident  in  the  transport 
profiles  (Fig.  5)  at  sections  A  and  B,  which 
are  near  tbe  Strait  of  Gibraltar,  and  show 
that  the  upper  portion  of  the  outflow  water 
begins  to  mix  with  the  overlying  North 
Atlantic  water  while  still  in  the  strait  (19). 

The  value  of  STwtd  then  decreases  rapidly 
to  about  36.6  psu  at  section  D,  which  is 
accompanied  by  a  significant  change  in  the 
transport  profiles:  from  sections  B  to  E, 
there  is  an  increase  of  about  0.6  X  106  m3 
s-1  in  the  westward  volume  transport  (out¬ 


flow)  of  water  in  the  density  range  ct6  =  27.3 
to  27.9  kg  m-3;  from  sections  E  to  B,  there 
is  a  complementary  decrease  in  the  eastward 
volume  transport  of  North  Atlantic  water  in 
the  density  range  ct9  =  27.0  to  27.2  kg  m-3. 
This  is  a  result  of  the  entrainment  of  North 
Atlantic  Central  water,  which  causes  the 
volume  transport  of  the  outflow  to  roughly 
double  and  the  value  of  STwtd  to  decrease 
rapidly  in  the  region  around  sections  C  and 
D.  Westward  of  section  E,  the  continued 
change  in  SXwtd  and  the  transport  profiles  is 
fairly  gradual,  with  STwtd  decreasing  to  about 
36.3  psu  at  section  H  in  the  extreme  western 
Gulf  of  Cadiz. 

The  rate  of  change  of  salinity  and  density 
along  the  path  of  the  outflow  is  strongly 
inhomogeneous.  Roughly  50%  of  the  net 
decrease  in  salinity  and  density  occurs  be¬ 
tween  sections  B  and  E,  which  is  only  about 
15%  of  the  total  path  length  in  the  Gulf  of 
Cadiz.  This  is  also  the  region  where  the 
outflow  begins  to  descend  the  continental 
slope  (Fig.  5A),  which  suggests  that  the  mix¬ 
ing  depends  on  bottom  topography  and,  by 
inference,  the  dynamics  of  the  current  itself. 

Elements  of  Outflow  Dynamics 

Laboratory  studies  provide  considerable  in¬ 
sight  into  the  dynamics  of  density-driven 
currents  (20,  21).  Even  though  no  single 
laboratory  experiment  includes  all  of  the 


Fig.  3.  Profiles  of  CTD-measured  salinity  and  XCP-measured  currents  along  the  outflow  axis  at 
sections  (from  right  to  left)  A,  C,  and  F,  The  current  profiles  are  viewed  as  if  looking  toward  the 
northeast  from  a  vantage  point  above  the  Gulf  of  Cadiz.  The  XCP  profiles  show  a  jetlike  current  of 
Mediterranean  outflow  water  with  a  typical  thickness  of  100  m  and  a  maximum  speed  from  sections 
A  to  D  of  -1  ms-1.  In  most  profiles,  the  current  maximum  was  near  the  middle  of  the  outflow  layer, 
and  there  were  generally  thick  shear  layers  in  the  interfacial  boundary  layer  between  the 
Mediterranean  water  and  the  inflowing  North  Atlantic  water,  and  in  the  bottom  boundary  layer.  The 
depth  of  the  outflow  increased  from  east  to  west  as  it  descended  the  continental  shelf  and  slope. 
The  outflow  current  changed  direction  by  about  90°  and  slowed  considerably  from  sections  C  to  F. 
The  width  of  the  outflow  layer  increased  from  about  10  km  at  section  A  to  about  45  km  at  section  F, 
so  that  the  volume  transport  within  the  outflow  layer  increased  by  about  a  factor  of  2  from  sections 
A  to  F.  The  salinity  profiles  show  a  substantial  decrease  in  the  outflow  salinity  going  from  sections 
A  to  F.  Within  the  outflow  layer,  the  salinity  (and  density)  profiles  typically  show  a  quasi-mixed  layer 
near  the  bottom  that  is  about  half  the  total  thickness  of  the  outflow  layer  and  a  highly  stratified 
interface  layer  of  comparable  thickness. 
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complexity  of  an  oceanic  outflow,  we  can 
nevertheless  deduce  some  of  the  important 
elements  of  outflow  dynamics  from  idealized 
examples.  A  simple  but  still  relevant  case  is 
that  of  a  source  of  constant  buoyancy  that 
feeds  a  density  current  that  is  descending  a 
straight  (nonrotating),  inclined  channel. 
The  current  accelerates  down  the  slope 
until  it  reaches  a  quasi-steady  state  in  which 
the  pressure  gradient  along  the  stream  is 
balanced  by  the  retarding  effect  of  bottom 
stress  and  entrainment  stress. 

The  importance  of  the  entrainment 
sttess  depends  on  the  internal  Froude  num¬ 
ber,  F  =  U/Vgfi  (21,  22),  where  U  is  the 
current  speed,  g  =  gAp/p0  is  the  reduced 
gravity,  g  is  the  acceleration  of  gravity,  Ap 
is  the  density  difference  between  the  densi¬ 
ty  current  and  the  ambient  water  (the 
density  anomaly),  p0  is  a  reference  density, 
and  h  is  the  thickness  of  the  density  cur¬ 
rent.  The  Froude  number  appears  in  many 
hydraulics  problems  and  has  several  inter¬ 
pretations,  including  that  it  is  the  ratio  of 
the  current  speed  to  the  phase  speed  of  long 
internal  waves.  If  F  <  1,  the  density  cur¬ 
rent  will  not  mix  significantly  with  the 
overlying  water.  In  that  case,  the  retarding 
stress  that  balances  the  pressure  gradient  is 
mainly  bottom  stress.  If  F  >  1,  the  entrain¬ 
ment  stress  may  provide  a  significant  part  of 
the  total  stress.  In  that  case,  the  density 
current  will  entrain  the  overlying  water  and 
consequently  lose  some  of  its  density  anom¬ 
aly.  This  suggests  that  F  is  a  crucial  internal 
flow  parameter  of  an  outflow.  Estimating  F 
in  a  straightforward  way  from  the  CTD  and 
XCP  data,  we  found  that  F  **  1  in  the  first 
50  to  100  km  of  the  outflow  path  and  is 
significantly  less  over  the  next  100  km. 


R#.  4.  A  profile  of  current  and  salinity  as  in  Fig. 
3  and  a  profile  of  dissipation  from  an  expand¬ 
able  dissipation  profiler  (XDP)  (right  panel). 
These  data  are  from  a  typical  station  along  the 
outflow  axis  at  section  C.  Dissipation  was  ex¬ 
tremely  high  in  the  saline  outflow  layer,  of  order 
10~*  W  nr3,  and  was  about  a  factor  of  1 04  less 
in  the  inflowing  Atlantic  water  above. 


To  understand  why  the  Froude  number  of 
the  outflow  might  exceed  1,  we  examined 
how  the  momentum  balance  changes  along 
the  path  in  response  to  variation  of  external 
parameters.  The  bottom  stress  on  a  turbulent 
current  is  often  approximated  by  a  quadratic 
drag  law  th  =  pQCD(J2,  where,  for  ocean 
bottoms,  CD  —  (1.5  to  3)  x  10-3,  depend¬ 
ing  on  the  roughness  of  the  terrain  and  the 
height  at  which  U  is  observed.  Under  this 
approximation,  it  is  easy  to  show  that  the 
Froude  number  of  a  quasi-steady  density 
current  will  exceed  1  if  dx  +  K  —  cD,  where 
dx  and  hx  are  the  along-stream  slopes  of  the 
bottom  depth  and  the  layer  thickness  (23). 
This  suggests  that  the  bottom  slope  and  the 
drag  coefficient  are  important  external  pa¬ 
rameters  of  a  density  current. 

The  terms  in  the  momentum  balance  can 
also  be  estimated  by  an  analysis  of  the  data. 
The  bottom  and  entrainment  stresses  can  be 
inferred  from  models  of  boundary  layer  flow. 
We  used  the  XCP  measurements  of  the 
current  in  the  lowest  5  to  10  m  of  the  bottom 
boundary  layer  to  estimate  the  bottom  stress, 
Th,  by  fitting  U(?)  to  the  logarithmic  form 
(U./k)  log(z/Zo),  where  U,  =  is  the 

friction  velocity,  which  is  to  be  determined, 
k  =  0.4  is  von  Karman’s  constant,  and  z0  is 
the  roughness  length  (24).  The  estimated 
bottom  stress  is  quite  large  and  shows  a 
systematic  variation  along  the  path  that  is 
roughly  proportional  to  the  variation  of  U2: 
Th  ~  1.5  Pa  at  sections  A  and  B,  has  a 
maximum  of  about  2.5  Pa  at  section  C,  and 
is  much  less,  about  0.3  Pa,  at  section  F. 

The  XDP  measurements  of  the  turbu¬ 
lence  intensity  were  used  to  estimate  the 
profile  of  turbulent  kinetic  energy  dissipa¬ 
tion,  e(r) .  This  too  can  be  used  to  infer 


stress,  with  the  assumption  of  a  simplified 
energy  budget,  e(z)  =  t(?)  U?,  where  t(z)  is 
the  turbulent  stress  and  Uz  is  the  vertical 
shear  of  the  horizontal  current,  which  can 
be  readily  estimated  from  XCP  data  (25). 
We  used  this  method  to  estimate  the  ratio 
of  the  turbulent  stress  in  the  interfacial 
layer,  Te,  to  the  bottom  boundary  layer 
stress,  T  =  TjTb,  and  found  that  the  mean 
value  over  29  samples,  distributed  mainly 
between  sections  A  and  D,  is  T  =  0.4  ± 
0.2,  where  the  uncertainty  is  twice  the 
standard  error.  Thus  the  bottom  stress  ap¬ 
pears  to  exceed  the  entrainment  stress  by 
about  a  factor  of  2  over  that  portion  of  the 
outflow  path.  There  appears  to  be  a  spatial 
variation  of  T  within  this  region,  although 
the  small  number  of  samples  together  with 
the  large  variance  of  T  render  this  somewhat 
tenuous.  The  mean  (and  median)  at  sections 
A  and  B  is  T  =  0.3  ±  0.1  (0.2)  (n  =  9); 
there  appears  to  be  a  maximum  at  section  C, 
where  T  =  0.9  ±  0.5  (0.6)  (n  =  11),  and  a 
small  decrease  at  section  D,  where  T  =  0.6 
±  0.6  (0.4)  (n  =  6).  The  number  of  XDP 
samples  from  sections  E  and  F  is  inadequate 
to  yield  a  stable  mean  of  T,  but  inspection  of 
the  dissipation  profiles  suggests  that  T  may 
decrease  further  to  the  west. 

The  pressure  gradient  averaged  over  the 
outflow  layer  can  be  written  as  Px  = 
g|(Ap/i)x  -I-  ApdJ,  where  d  is  the  bottom 
depth  and  the  derivative,  indicated  by  the 
subscript  x,  is  taken  along  the  direction  of 
the  outflow.  The  first  term  represents  the 
internal  (baroclinic)  pressure  gradient, 
which  arises  from  variations  in  density  or 
layer  thickness  along  the  outflow  path.  This 
term  is  dominant  within  the  Strait  of 
Gibraltar  proper  and  drives  the  two-layer 


Flfl.  5.  (A)  Volume  transport  across 
each  of  the  sections  shown  as  a 
function  of  depth.  The  outflow  wa¬ 
ter  is  blackened.  (B)  Volume  trans¬ 
port  shown  as  a  function  of  density 
(we  use  density  rather  than  salinity 
because  density  is  monotonically 
increasing  with  depth,  while  salini¬ 
ty  has  a  local  minimum  in  the  At¬ 
lantic  water  just  above  the  outflow 
layer).  Note  the  large  change  in  the 
transport  profile  as  a  function  of 
density  between  sections  B  and  E. 
This  is  evidence  of  entrainment  by 
the  outflow  of  North  Atlantic  Central 
water  from  the  density  range  <x„  = 
27.0  to  27.2  kg  ITT3. 
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exchange  flow.  Mediterranean  water  is 
forced  westward  through  the  strait  at  a 
speed  of  about  1ms-1  [although  the  speed 
is  highly  dependent  on  the  tidal  phase 
within  the  strait;  see  (19)  for  a  detailed 
discussion  of  the  complex  time-dependent 
hydraulic  control  phenomenon  within  the 
strait).  After  the  main  sill  at  Camarinal,  the 
Mediterranean  water  overflows  another  mi¬ 
nor  sill  at  6°10'W,  which  is  very  close  to 
our  section  A,  and  begins  a  monotonic 
descent  into  the  Gulf  of  Cadiz. 

Mediterranean  Outflow  Dynamics 

Once  the  outflow  enters  the  Gulf  of  Cadiz, 
the  pressure  gradient  is  increasingly  domi¬ 
nated  by  the  second  term,  the  topographic 
term,  and  the  outflow  is  very  strongly  af¬ 
fected  by  bottom  topography,  as  earlier 
investigators  have  emphasized  in  different 
ways  (8,  10-13).  Our  description  of  the 
outflow  dynamics  is  divided  into  three  seg¬ 
ments  that  are  defined  by  the  relation  of  the 
outflow  current  to  the  bottom  topography. 

Channel  flow.  At  sections  A  and  B,  the 
Mediterranean  outflow  is  constrained  with¬ 
in  a  narrow  channel  ~  10  km  wide  that  runs 
west-southwest.  This  channel  appears  to  be 
an  extension  of  the  Strait  of  Gibraltar  that 
has  been  cut  out  of  the  continental  shelf  by 
the  outflow  current.  The  along-stream  slope 
of  the  channel  bottom  is  typical  of  a  con¬ 
tinental  shelf,  dx  ~  (1  ±  0.5)  x  UT\ 
where  the  large  uncertainty  indicates  the 
variation  on  the  scale  of  the  outflow  width. 
The  rate  of  change  of  the  layer  thickness  is 
similar;  thus,  dx  +  hx  S  CD.  The  observed 
Froude  number  is  also  Si,  and,  as  we  have 
seen,  there  is  some  mixing  in  this  region, 
although  it  is  not  particularly  strong. 

The  pressure  force  per  unit  area  along 
the  direction  of  the  outflow  is  estimated  to 
be  hPx  =  hgAp(dx  +  hx)  *  3  ±  1.5  Pa, 
where  Ap  =  1.5  kg  m-3.  This  is  a  large 
force  that  would,  if  unbalanced  by  a  stress, 
cause  the  outflow  current  to  accelerate  from 
about  1  ms'1  at  section  A  to  about  1.4  m 
s-1  by  section  B.  No  such  acceleration  is 
observed,  evidently  because  this  pressure 
force  is  largely  balanced  by  the  sum  of  the 
bottom  stress,  which  is  rh  *  1.5  Pa  in  this 
region,  estimated  from  the  fit  to  XCP  data, 
and  the  entrainment  stress,  which  is 
0.5  Pa,  estimated  from  the  ratio  t Jtb  that 
was  inferred  from  XDP  and  XCP  data.  The 
combined  bottom  and  entrainment  stress 
gives  a  retarding  stress  of  t  *  2  Pa,  which 
appears  sufficient  to  roughly  balance  the 
pressure  force. 

This  first  region  of  the  outflow  can  thus  be 
characterized  as  a  more  or  less  steady  channel 
flow  that  has  a  Froude  number  of  —1.  The 
bottom  stress  provides  most  of  the  retarding 
stress,  evidently  because  the  observed  slopes 
do  not  exceed  the  drag  coefficient. 


Initial  descent  of  the  continental  slope.  The 
outflow  reaches  the  edge  of  the  continental 
shelf  just  west  of  section  B.  It  flows  directly 
across  the  much  steeper  continental  slope 
near  section  C  where  dx  =  (11  ±  2)  X 
10~3,  which  is  considerably  larger  than  the 
drag  coefficient.  The  observed  internal 
Froude  number  reaches  a  maximum  of 
about  1.3,  and  entrainment  is  much  stron¬ 
ger  here  than  just  upstream. 

The  momentum  balance  is  also  quite 
different.  The  steep  bottom  slope  causes  a 
very  large  pressure  force  where  the  outflow 
first  crosses  onto  the  continental  slope: 
hgApdx  ~  13  Pa,  where  Ap  =  1.4  kg  m-3  is 
the  initial  value  of  the  density  anomaly. 
The  current  accelerates  to  a  maximum 
speed  of  about  1.3  m  s'1  at  section  C, 
which  is,  again,  very  much  less  than  what 
would  be  expected  from  this  pressure  force. 
There  appear  to  be  three  processes  that 
limit  the  downhill  acceleration  of  the  cur¬ 
rent  in  this  crucial  region. 

1)  Turbulent  stresses  that  arise  from  en¬ 
trainment  and  bottom  friction  are  very 
large  and  produce  intense  dissipation 
throughout  the  entire  thickness  of  the  out¬ 
flow  (Fig.  4)  •  Many  individual  XDP  profiles 
show  two  distinct  dissipation  maxima:  an 
upper  maximum  in  the  interfacial  boundary 
layer  between  the  outflow  and  the  overlying 
Atlantic  water  and  a  lower  maximum  in  the 
bottom  boundary  layer.  The  total  stress  that 
was  estimated  from  the  combined  XCP  and 
XDP  data  is  about  5  Pa  at  section  C,  with 
roughly  equal  contributions  from  entrain¬ 
ment  stress  and  bottom  stress.  The  com¬ 
bined  stress  represents  a  very  significant 
brake  on  the  outflow  current,  but  it  is  only 
about  one-third  of  the  initial  pressure  force. 

2)  Entrainment  has  an  additional,  and 
in  some  respects  decisive,  effect  on  the 
outflow  by  reducing  its  density,  and  thus, 
the  density  anomaly.  The  outflow  density  is 
reduced  to  about  cre  =  27.75  kg  m-5  by 
section  D,  and  the  density  anomaly  is  re¬ 
duced  to  about  Ap  «=  0.7  kg  m-3,  which  is 
a  little  less  than  half  of  its  value  at  section 
B.  This  halves  the  pressure  force  on  the 
outflow. 

3)  Finally,  the  earth’s  rotation  also  lim¬ 
its  the  downhill  acceleration  of  the  outflow: 
the  Coriolis  force  turns  the  current  to  the 
right  (looking  in  the  direction  of  the  flow). 
Instead  of  flowing  directly  down  the  steep¬ 
est  slope,  the  outflow  turns  and  flows  more 
nearly  parallel  to  isobaths  (Fig.  1C).  By 
following  a  column  that  was  moving  from 
section  B  to  D,  we  estimated  that  the 
current  turns  at  —90°  per  7  hours,  which  is 
comparable  to  the  turning  rate  of  an  inertial 
current  at  this  latitude  (the  Coriolis  param¬ 
eter/  =  90°  per  5.1  hours)  (26). 

The  initial  descent  of  the  continental 
slope  thus  has  the  character  of  a  geostrophic 
adjustment  process  in  which  the  outflow 


current  responds  to  an  increased  bottom 
slope  (27).  The  increased  pressure  force 
accelerates  the  current  enough  to  raise  the 
Froude  number  above  1,  and  entrainment 
becomes  intense.  The  consequent  loss  of 
density  anomaly,  together  with  the  effect  of 
geostrophic  adjustment,  limits  this  region  of 
strong  mixing  to  only  a  few  tens  of  kilome¬ 
ters  or,  following  the  current,  only  a  few 
hours.  However,  in  this  brief  interval,  the 
density  of  the  outflow  is  reduced  enough 
that,  even  if  no  further  mixing  occurred,  the 
Mediterranean  outflow  would  settle  into  the 
thermocline  rather  than  the  deep  sea. 

Damped  geostrophic  flow.  As  the  outflow 
reaches  sections  E  and  F,  it  turns  to  flow 
nearly  parallel  to  the  bottom  topography 
(Fig.  1C).  The  current  speed  is  greatly 
reduced  to  about  U  =  0.4  ms-1,  and  the 
current  is  then  in  a  near  geostrophic  bal¬ 
ance  in  which  the  Coriolis  acceleration  is 
roughly  balanced  by  the  across-stream  topo¬ 
graphic  acceleration:  fU  ~  gdy,  where  Ap  = 
0.5  kg  m-5  and  dy  =  (8  ±  3)  X  10“  *  is  the 
bottom  slope  normal  to  the  current.  If  the 
Froude  number  of  a  geostrophic  outflow  is 
to  exceed  1 ,  then  for  the  observed  thickness 
and  density  anomaly,  the  bottom  slope 
would  have  to  satisfy  dy  >  /  \ZK/£  ~  13  x 
10“ 3,  which  it  does  not  near  section  F  (28). 
The  maximum  observed  Froude  number  is 
about  0.3,  well  below  the  level  at  which 
significant  entrainment  might  be  expected. 
The  comparatively  tranquil  character  of  the 
flow  is  evident  in  the  XCP  and  XDP  profiles 
(Figs.  3  and  4),  which  show  that  the  inter- 
facial  boundary  layer  has  a  much  reduced 
thickness  and  turbulence  intensity.  This  sug¬ 
gests  that  mixing  between  the  outflow  and 
Atlantic  water  is  considerably  reduced  at 
sections  E  and  F,  which  is  consistent  with 
the  much  slower  rate  of  decrease  of  STwld 
along  the  path  in  this  region  (Fig.  2C). 

The  turbulence  in  the  bottom  boundary 
layer  is  also  much  less  intense  in  this  re¬ 
gion,  and  the  bottom  stress  is  about  0.3  Pa 
as  estimated  from  XCP  profile  data  along 
section  F.  This  weaker  but  still  significant 
bottom  stress  may  play  a  crucial  role  in  the 
longer  term  evolution  of  the  outflow  by 
providing  the  damping  (dissipation)  that 
allows  the  otherwise  geostrophic  current  to 
descend  the  continental  slope  without  ac¬ 
celerating.  The  effect  of  bottom  friction  on 
a  nearly  geostrophic  current  can  be  inferred 
from  the  ratio  of  the  bottom  stress  to  the 
across-stream  (geostrophic)  topographic  ac¬ 
celeration,  which  is  known  as  the  Ekman 
number,  K  =  (rh/p h)/(gdy),  where  h  =  100 
m  and  Ap  and  d,  are  as  before.  Thus,  K  « 
0.2,  which  is  the  angle  of  the  current,  in 
radians,  with  respect  to  the  isobaths.  As  the 
outflow  moves  from  section  F  to  Cape  Saint 
Vincent,  a  distance  L  =  200  km,  it  can 
descend  an  additional  depth  LKdy  =  200  m 
without  accelerating. 
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Our  hydrographic  data  show  that  the 
outflow  continues  westward  along  the  con¬ 
tinental  slope  as  a  damped  geostrophic  cur¬ 
rent  that  slowly  crosses  isobaths.  There  is 
additional  mixing,  although  the  rate  is 
much  less  than  that  found  around  section 
C,  and  when  the  outflow  reaches  Cape 
Saint  Vincent  and  the  open  North  Atlan¬ 
tic,  it  has  properties  that  make  it  neutrally 
buoyant  at  depths  of  about  800  to  1300  m: 
T  =  12°  ±  1°C,  S  =  36.4  ±  0.1  psu,  and 
ct9  =  27-6  ±  0.1  kg  m-5. 

Other  North  Atlantic  Outflows 

The  lower,  denser  portion  of  the  outflow 
floats  off  of  the  bottom  in  the  western  Gulf 
of  Cadiz  and  forms  the  saline  tongue  that 
extends  west  across  the  North  Atlantic 
basin  (Fig.  1A)  (29).  The  upper  portion  of 
the  outflow  continues  north  for  some  dis¬ 
tance  along  the  Iberian  Peninsula  and 
merges  with  a  poleward  eastern  boundary 
undercurrent  that  carries  the  saline  Medi¬ 
terranean  water  toward  the  inflow  of  the 
Norwegian  Sea,  just  west  of  Ireland.  The 
relatively  high  salinity  of  the  water  that 
flows  into  the  Norwegian  Sea  must  enhance 
the  density  of  the  very  cold  waters  that  are 
produced  there  during  the  winter.  Cold, 
dense  water  from  the  Norwegian  Sea  forms 
two  important  outflows  (30-32)  that  form 
the  deep  and  bottom  water  in  the  North 
Atlantic.  By  contributing  to  the  salinity  of 
the  Norwegian  Sea  inflow  (33),  the  warm, 
dry  climate  over  the  Mediterranean  Sea 
may  enhance  the  cold,  thermally  dominat¬ 
ed  circulation  of  the  deep  ocean  (14). 

The  qualitative  pattern  of  |he  mixing 
and  dynamics  found  in  the  Mediterranean 
outflow — intense  mixing  during  the  initial 
descent  of  the  continental  slope  followed  by 
a  long  stretch  of  nearly  geostrophic  flow- 
may  be  discerned  as  well  in  historical  ob¬ 
servations  of  the  Norwegian  Sea  outflow 
through  the  Denmark  Strait  (32).  This 
outflow  begins  at  the  sill  west  iof  Iceland 
with  properties  that  make  it  less  |ense  than 
pure  Mediterranean  outflow  water:  T  = 
0°C,  S  =  34.9 psu,  and <re  =  28.05  kg m~3. 
It  flows  southward  along  the  continental 
slope  of  Greenland  and  entrains  warmer, 
presumably  overlying,  water  somewhere  in 
the  northern  Irminger  Sea,  which  is  evi¬ 
dent  because  it  warms  1°  to  2°C.  This 
outflow  loses  only  about  0.1  kg  m“3  of  its 
density  by  entrainment  and  remains  just 
dense  enough  to  become  bottom  water  in 
the  North  Atlantic.  The  Denmark  Strait 
outflow  continues  south  for  thousands  of 
kilometers  along  the  continental  slope  of 
North  America  as  the  deepest  and  densest 
portion  of  a  geostrophic  current  called  the 
Deep  Western  Boundary  Current.  This  cur¬ 
rent  mixes  only  very  slowly  with  its  sur¬ 
roundings  (34)  and  retains  its  chemical 


identity  at  least  as  far  south  as  the  tropics. 

Models  of  the  general  circulation  that 
simulate  a  deep  ocean  that  changes  over 
time  will  probably  need  to  include  some 
representation  of  time-varying  marginal  sea 
outflows.  There  are  at  least  two  somewhat 
different  problems  posed  by  this:  first,  to 
specify  or  simulate  the  seawater  properties 
of  the  marginal  seas,  and  second,  to  allow 
the  marginal  seawater  to  flow  out  into  the 
open  ocean  in  a  realistic  manner.  Our  study 
gives  some  indications  of  what  may  be 
required  for  this  second  part.  The  geo¬ 
strophic  portion  of  an  outflow  current  and 
the  accompanying  gradual  changes  in  out¬ 
flow  properties  would  seem  to  present  no 
special  demands.  However,  the  rapidly 
varying  portion  of  an  outflow  path,  for 
example,  the  first  100  km  of  the  Mediter¬ 
ranean  outflow  where  the  transport  dou¬ 
bles,  will  pose  resolution  and  parameteriza¬ 
tion  problems  that  appear  to  be  unique  to 
outflows.  These  particular  portions  of  the 
outflows,  which  we  think  are  localized 
around  steep  topography  near  the  marginal 
sea,  may  have  to  be  treated  by  a  submodel 
(35)  that  has  enhanced  vertical  resolution 
and  mixing  parameterizations  which  are 
appropriate  to  density  currents. 
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